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PREPACE 

This  Mtmorandum  Is  concerned  with  the  application 
of  Invariant  Isibeddlng  to  the  study  of  rarefied  gas 
dynamics.  A  new  formulation  of  the  problem  of  linearised 
Couette  flow  Is  presented  using  Ideas  which  were  first 
developed  In  astrophysics  and  neutronrtransport  ^eory. 
This  study  tiiould  be  of  Interest  to  qpeoiallata  in  the 
fields  of  aerodynamics  and  heat  transfer  who  are 
concerned  with  the  calculation  of  flow  fields  In  the 
upper  atmosphere*  as  well  as  to  astrophysicists  and 
nuclear  i^slclsts. 


■SUMMAHy 


This  Namorandtaa  applies  the  teehnlqpes  of 
Invariant  Inbeddlng  to  the  study  of  rarefied  gas  flows. 
The  problem  of  linearised  Couette  flow  Is  Investigated, 
and  It  Is  shown  how  the  aasusptlon  of  the  Krook 
scattering  model  results  In  a  foraailatlon  idileh  Is 
similar  to  that  obtained  In  radiative  transfer  for 
cmiservatlve  Isotropic  scattering  In  a  plane-parallel 
atmosphere. 

By  a  slxople  emaneratlon  of  i^slcal  processes, 
the  nonlinear  Integral-dlffexentlal  eqoatlcn  governing 
the  reflectlmi  function  la  obtained,  and  a  suitable 
transformatlm  Is  shown  to  render  this  equation  amenable 
to  numerical  ccnputatlcm. 
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1.  INTRODUCTION 

The  use  of  Invazlance  principles  and  Invariant 
Imbedding  techniques*  originating  in  the  work  of 
Ambarziimian  and  Chandrasekhar*  has  provided  new  insights 
and  conqputatlonal  solutions  to  various  problems  in  the 
fields  of  radiative  transfer  and  neutron  transport 
theory  [ iJ .  Let  us  indicate  how  these  same  techniques 
can  be  applied  to  representative  problons  in  the  kinetic 
theory  of  gases,  nie  possibility  of  this  is  suggested 
by  the  essential  similarity  of  all  particle  processes: 
the  b^avior  of  neutrons  or  photons*  or  the  structureless 
spherically  symmetric  molecule  of  the  classical  kinetic 
theory.  For  the  present*  we  shall  restrict  our  area  of 
interest  to  linearised  scattering  models*  finite 
collision  frequencies*  and  relatively  simple  gecsoetries. 
We  leave  to  future  study  the  application  of  invariance 
techniques  to  nonlinear  scattering  isodels  and  to  eases 
in  which  the  Boltzmann  gain  and  loss  operators  are  rnily 
conditionally  convergent.  A  classic  problem*  that  of 
plane  shear  flow  between  two  Infinite  flat  plates*  will 
be  used  to  present  the  ideas  of  invariant  labeddlag  in 
the  context  of  rarefied  gas  dynasdes. 
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2.  KINEriC-THBORY  APPROACH 

Various  klnetlc~theory  approaches  to  the  problem  of 
plane  shear  flow  are  described  In  some  detail  by  Villlls 
[2].  The  physical  picture  Is  the  following. 

An  Infinite  plane  wall  moves  with  a  small  constant 
velocity  In  Its  plane  at  a  distance  d  from  a  fixed 
wall.  Both  walls  are  maintained  at  the  same  temperature* 
and  may  consist  of  either  similar  or  dissimilar  materials. 

It  Is  desired  to  find  the  gas  velocity  and  the  stx^ss  at 
the  moving  wall  as  a  function  of  the  degree  of  rarefaction* 
which  is  measured  by  an  appropriate  Knudsen  nundt>er.  In 
the  body  of  the  gas*  we  shall  use  the  single  relaxation 
time  model  of  Krook  [3^  in  Its  linear  version.  Linearisation 
Is  made  possible  by  requiring  that  the  Mach  nunflser  of  the 
moving  wall  be  much  less  than  unity*  and  then  using 
standard  perturbation  technlq[ues. 

We  define  the  following  quantities: 

(2.1)  f 

u*v*w  ■  VE  UiVE  VE  1* 

The  nenllnsar  Ksoek  (3I  •qpwtlen  for  this  ono- 
tflmantioiial  goomotxy  can  be  written 


-3- 

(2.2)  ^  exp(-  h(c  -  V)^)] . 

The  linearized  version  of  this  equation*  obtained  by 
setting  V„  «  1,*  (where  V  Is  the  velocity  of  the 
moving  wall)  Is 

(2.3)  ^  ■  -  6rJf  -  n(y)^/^  exp(-  hc®)(l  +  2hv  V)] . 

If  f  -  nQ(  ^)3/2  «<p(_  hoC=)  +  f  where 

fl  ■  ®xp(—  hQC^)v^(u*x),  then  the  equation  for 

the  perturbation  distribution  function  may  be  written 

udf,  r  r  r+oo  p  “I 

(*•")  +  JfjYdudvdi.. 

i-  * 

The  boxmdaxy  conditions  corresponding  to  purely  diffuse 
reflection  from  the  wall  are 

(2.5)  (a)  Fixed  plate:  ft(x  «  0)  «  0* 

(b)  Moving  plate:  f^(x  ■  d)  «  e"®^»2  v  M,^ 

where  M,,  «  VF  V^. 

It  Is  possible  to  develop  sn  equation  for  tiie  function 
which  Is  closely  related  to  the  radiative  transport 
equation  for  Isotropic  scattering  In  a  plane  slab*  l.e.« 

(*.6)  J"** 

The  diffuse  boundary  eondltlons  aseeolated  «ltb  bbls 
ev*stlen  are 

(2.T)  y  -  0,  -  0,  y  -  1,  •*"  -  a\lK 

This  Is  equivalent  to  xufiiiring  that  the  Nash 
nisriber  be  everywhere  aaall.  m  Is  proportional  to  the 
Nah  nunfeer  of  the  wall.  ^ 


Willis  has  obtained  a  nianerlcal  solution  to  equation  (2.6) 
by  recasting  It  In  a  form  similar  to  the  Milne  Integral 
equation  of  radiative  transfer.  In  addition,  he  has 
shown  that  stralghtforvmrd  Neumann  Iteration  Is  suitable 
only  for  the  very  rarefied  case,  corresponding  to  t  <  1. 
Since  each  subsequent  calculation  in  the  Neumann  iteration 
scheme  corresponds  to  higher-order  collisions.  It  Is  not 
surprising  that  It  is  rather  slowly  convergent  as  the 
collision  frequency  is  Increased.  Dils  Is  quite  similar 
to  the  state  of  affairs  In  radiative  transfer,  where  It 
is  known  that  a  calculation  based  on  successive  absorption 
and  scattering  processes  Is  ill  suited  for  c(wnputatlon  as 
slab  thickness  Is  increased,  or  conversely,  as  the  optical 
mean  free  path  is  decreased. 

3.  INVABIANT  DBEDDINQ  APPROACH 

Ohe  invariant-imbedding  approach  permits  us  to 
cmoentrate  our  attention  on  ^e  valuer  of  quantities  at 
the  boundaries.  This  is  particularly  relevant  to  problcsui 
in  heat  transfer  and  aerodynamlosa  where  values  of 
certain  molecular  fluxes  at  the  boundary  are  related  to 
such  maorosoopio  observables  as  energy  transport  and 
momentum  transport  to  the  walls.  Once  values  have  been 
obtained  for  the  distribution  funotion  at  the  bountery# 
however^  it  is  possible  to  oaloulete  values  of  the  distri¬ 
bution  funotion  within  the  body  the  gas  by  relatively 
simple  teohniqtws. 


The  Introduction  of  other  bovindary  conditions 
besides  the  presently  chosen  diffuse  re-emlsslon  Is 


postponed  so  that  we  may  focus  on  the  essentials  of  the 
technique.  In  addition^  It  should  be  noted  that  our 
choice  of  one  moving  wall  and  one  fixed  wall  Is  someidiat 
different  from  the  usual  Couette  flow  problem^  where  the 
two  walls  move  In  opposite  directions.  This  choice  was 
made  because  It  obviated  the  need  for  the  calculation  of 
transmission  functions.  ‘Hie  method  of  linearization  and 
the  diffuse  boundary  conditions  at  the  fixed  wall  result 
In  the  fixed  wall  acting  as  a  "sink"  for  perturbing 
molecules.  In  the  sense  that  It  anlts  only  particles 
vdilch  possess  the  equlllbrlvun  distribution  function. 

Consider  Fig.  1.  A  x*ef lection  function 
p(d,u,u^,v,v^,w,w^)  Is  defined  vAilch  relates  Incoming 
moleciiles  at  d  to  outgoing  molecules  at  d: 

.  r  +00  r  +00  r  +00 

(3.1)  •  J  J  J  p(d,u,u^,v,v^,w,w^) 

0  -00  -00 


u»o 


Fig.  1 


Linearity  of  the  transport  equation  pe’Tnits  us  to 
consider  a  Dirac  6—functlon  Incident  distribution 
centered  about  Uq,Vq,Wq.  Thus  p(d,u,UQ,v,VQ,w,WQ)  is 
the  outgoing  distribution  function  of  particles  at  d 
due  to  an  incident  monochromatic  beam  at  d  vrLth  velocity 
Uo^VojWq.  We  observe  that  the  sink-like  character  of  the 
fixed  diffuse  wall  guarantees  that  all  outgoing  perturbing 
pa]:*ticles  at  the  moving  vail  are  due  solely  to  the 
disturbance  generated  by  it.  Simple  particle-co\mtlng 
techniques  are  utilized  to  relate  the  reflection  function 
p(d,u,UQ,v,VQ,w,WQ)  at  d  to  the  reflection  function 
p(d  +  A,u,Uq,v,  Vq»w,Wq)  at  d  +  A,  We  recognize  that 
the  original  disturbing  beam  of  molecxiles  at  d  +  A  is 
modified  by  Interactions  within  A,  and  the  resulting 
incoming  distribution  function  at  d  undergoes  the  same 
sort  of  processes  at  d  as  the  original  beam  underwent 
at  d  +  A.  A  new  problem  is  initiated  at  d,  which 
differs  from  the  problem  initiated  at  d  +  A  by  the  fact 
that  the  incident  distribution  is  continuous  rather  than 
moxiochroeiatie . 

we  introduce  the  mean  free  path  aa  the  vanXt  of 
length,  and  use  t  >  5  d  rather  than  d  for  the 

q>aoin8  betwemi  plates. 

Ihe  apirit  of  the  teohnlqfue  la  exhibited  In  the 
alai»le  partlole-oountlng  approach.  Vo  taxas  of  order 
(A)*,  p(t  t  A,u,Uq,v,Vq,w,Wq)  la  e«ual  to  the  roaolta 


of  the  following  processes: 

(a)  The  Incoming  monochromatic  beam  is  diminished 
by  collisions  in  A,  reflected  at  f,  and  then 
diminished  again  by  collisions  in  A. 

(b)  Ihe  incoming  monochromatic  beam  is  scattered 
within  A  from  Uq,Vq,Wq  to  u,v,w. 

(c)  ^e  incoming  beam  is  reflected  at  r,  and 
then  scattered  into  u,v,w  within  A. 

(d)  The  incoming  beam  is  scattered  within  A, 
then  reflected  at  r  into  u,v,w. 

(e)  Hie  incoming  beam  is  reflected  at  r, 
scattered  in  A,  and  then  reflected  again  at  t. 

Taking  account  of  these  interactions,  we  have 

(3.2)  p(T+d,U,UQ,V,VQ,W,WQ) 

-  [1  -  ^]p{t,u,Uq,v,Vq,w,Wq)[1  -  ^1  +  ■  •'^VqA 

u  V  V  V  uq  U 

f  ^p(‘^»w,UQ,t,VQ,w,WQ)tdad^  ♦ 

(▼)  ^  0  -00  -<0 

cOTiseq^ence  of  the  llxiearlsed  Krook  scattering 
model  Is  that  ATf2(u,v,w)/u  Is  the  number  of  pertloles 
idilch  ere  soattereo  out  of  the  element  veleoi^  space 
centered  about  u,v,w,  and 

**  f  Vi*^'**i*'i 

is  the  nMber  of  particles  scattered  into  this  elemsnt# 
all  in  an  infinitesimal  thickness  A. 
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y'  00 j'  +00  j  +00 


(ir) 


37? 


-00  -00 


2  V  V* 

e“^'  —  p(T,u,ti,v,^,w,<ir)(l4d^dilr 

u 


+  7  \372  ^  j  j  J  p(T,u,u,v,v,w,w)e“^  vdudvdS? 
\  **  / 


CO  J  OOJ  CO 

0  -00  -00 
00  r  +00  r  +oo 


/CO  r  +00  p  +00 

J  j  p(T,u,UQ,v,VQ,w,WQ)vdudvdw. 


0  -00  -oo 


a!he  reqtilranents  of  Sec.  2  suggest  the  transformation 

o  -52 

(3.3)  p(‘r*u,UQ,v,VQ,w,WQ)  - e”^  S(t,u,Uq), 


Letting  A  -*  0«  we  obtain  the  following  symmetric  fom 
for  the  one-dimensional  reflection  function  S(d«u,UQ): 


(3.4)  +  (i  +  i)S(T,u,UQ) 


S(t,U.U(,)  -^  +  S(T,U,Ug)^ 

i/"  S(t,u.»5) 

0  0  ^0 


The  initial  condition  associated  with  (3*4)  is 
S(0,u,Uq)  «  0«  resulting  from  the  observation  that 
t  ■  0  corresponds  to  the  free  molecule  limit  when 
there  are  no  collislwis  within  d  to  reflect  molecules 
from  Uq  Into  u. 

(Die  last  equation  closely  resmSbles  the  radiative 
transport  equation  for  Isotropic  soattarlng  In  a  plana 
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slab  that  Chandrasekhar  also  obtains  by  invariance 
principles  [4].  A  significant  difference^  however, 
between  radiative  transport  and  kinetic  theory  is  that 
in  kinetic  theory,  incoming  distribution  functions  are 
not  generally  of  a  delta-function  character,  so  that  it 
is  necessary  to  Integrate  over  all  Incoming  velocities 
in  order  to  obtain  the  distribution  function  at  the  wall. 

The  conQ)utational  solution  of  (3.4)  has  been 
readily  obtained  using  techniques  which  were  originally 
developed  for  problems  in  radiative  transfer.  In 
forthcoming  publications  we  shall  discuss  this,  consider 
the  inclusion  of  nondlffuse  particle-surface  interaction 
at  the  moving  wall,  and  present  asymptotic  observations 
which  are  relevant  to  the  problems  of  neaa^-f  ree  molecule 
flow  and  slip  flow. 
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